function [f,fx11,fx11p,fy,fyp,fx12] = getModelDeriv_EP(nx,mx,myx,ny,params,N,h0,g0,x,y,t,getDeriv)
 
%Unfold params
BETTA= params.BETTA;
B= params.B;
CHI= params.CHI;
CHI0= params.CHI0;
THETA= params.THETA;
DELTA= params.DELTA;
ALFA= params.ALFA;
PHI= params.PHI;
PHIzero= params.PHIzero;
ZI= params.ZI;
KAPAw= params.KAPAw;
ETA= params.ETA;
RHOR= params.RHOR;
PHIpai= params.PHIpai;
PHIpai_1= params.PHIpai_1;
PHIy= params.PHIy;
PHIy_1= params.PHIy_1;
PHIc= params.PHIc;
PHIc_1= params.PHIc_1;
PHIl= params.PHIl;
NU= params.NU;
U0= params.U0;
U0d= params.U0d;
OMEGAz= params.OMEGAz;
OMEGAd= params.OMEGAd;
OMEGAn= params.OMEGAn;
OMEGAp= params.OMEGAp;
OMEGAa= params.OMEGAa;
RHOz= params.RHOz;
RHOd= params.RHOd;
RHOn= params.RHOn;
RHOp= params.RHOp;
RHOa= params.RHOa;
Kss= params.Kss;
OUTPUTss= params.OUTPUTss;
PAIss= params.PAIss;
MUZss= params.MUZss;
Css= params.Css;
AA= params.AA;
lss= params.lss;
Rss= params.Rss;
Wss= params.Wss;
 
% Current values in the model;
% x(N,mx) y(N,ny);
% x;
c_ba1 = x(t,1)+h0(1,1);
muz_cu = x(t,2)+h0(2,1);
d_cu = x(t,3)+h0(3,1);
n_cu = x(t,4)+h0(4,1);
paistar_cu = x(t,5)+h0(5,1);
a_cu = x(t,6)+h0(6,1);
% xp;
c_ba1p = x(t+1,1)+h0(1,1);
muz_cup = x(t+1,2)+h0(2,1);
d_cup = x(t+1,3)+h0(3,1);
n_cup = x(t+1,4)+h0(4,1);
paistar_cup = x(t+1,5)+h0(5,1);
a_cup = x(t+1,6)+h0(6,1);
 
% y
c_cu = y(t,1)+g0(1,1);
r_cu = y(t,2)+g0(2,1);
pai_cu = y(t,3)+g0(3,1);
w_cu = y(t,4)+g0(4,1);
l_cu = y(t,5)+g0(5,1);
output_cu = y(t,6)+g0(6,1);
mc_cu = y(t,7)+g0(7,1);
evf_cu = y(t,8)+g0(8,1);
vf_cu = y(t,9)+g0(9,1);
p1_cu = y(t,10)+g0(10,1);
p2_cu = y(t,11)+g0(11,1);
p3_cu = y(t,12)+g0(12,1);
p4_cu = y(t,13)+g0(13,1);
p5_cu = y(t,14)+g0(14,1);
p6_cu = y(t,15)+g0(15,1);
p7_cu = y(t,16)+g0(16,1);
p8_cu = y(t,17)+g0(17,1);
p9_cu = y(t,18)+g0(18,1);
p10_cu = y(t,19)+g0(19,1);
p11_cu = y(t,20)+g0(20,1);
p12_cu = y(t,21)+g0(21,1);
p13_cu = y(t,22)+g0(22,1);
p14_cu = y(t,23)+g0(23,1);
p15_cu = y(t,24)+g0(24,1);
p16_cu = y(t,25)+g0(25,1);
p17_cu = y(t,26)+g0(26,1);
p18_cu = y(t,27)+g0(27,1);
p19_cu = y(t,28)+g0(28,1);
p20_cu = y(t,29)+g0(29,1);
p21_cu = y(t,30)+g0(30,1);
p22_cu = y(t,31)+g0(31,1);
p23_cu = y(t,32)+g0(32,1);
p24_cu = y(t,33)+g0(33,1);
p25_cu = y(t,34)+g0(34,1);
p26_cu = y(t,35)+g0(35,1);
p27_cu = y(t,36)+g0(36,1);
p28_cu = y(t,37)+g0(37,1);
p29_cu = y(t,38)+g0(38,1);
p30_cu = y(t,39)+g0(39,1);
p31_cu = y(t,40)+g0(40,1);
p32_cu = y(t,41)+g0(41,1);
p33_cu = y(t,42)+g0(42,1);
p34_cu = y(t,43)+g0(43,1);
p35_cu = y(t,44)+g0(44,1);
p36_cu = y(t,45)+g0(45,1);
p37_cu = y(t,46)+g0(46,1);
p38_cu = y(t,47)+g0(47,1);
p39_cu = y(t,48)+g0(48,1);
p40_cu = y(t,49)+g0(49,1);
% yp;
c_cup = y(t+1,1)+g0(1,1);
r_cup = y(t+1,2)+g0(2,1);
pai_cup = y(t+1,3)+g0(3,1);
w_cup = y(t+1,4)+g0(4,1);
l_cup = y(t+1,5)+g0(5,1);
output_cup = y(t+1,6)+g0(6,1);
mc_cup = y(t+1,7)+g0(7,1);
evf_cup = y(t+1,8)+g0(8,1);
vf_cup = y(t+1,9)+g0(9,1);
p1_cup = y(t+1,10)+g0(10,1);
p2_cup = y(t+1,11)+g0(11,1);
p3_cup = y(t+1,12)+g0(12,1);
p4_cup = y(t+1,13)+g0(13,1);
p5_cup = y(t+1,14)+g0(14,1);
p6_cup = y(t+1,15)+g0(15,1);
p7_cup = y(t+1,16)+g0(16,1);
p8_cup = y(t+1,17)+g0(17,1);
p9_cup = y(t+1,18)+g0(18,1);
p10_cup = y(t+1,19)+g0(19,1);
p11_cup = y(t+1,20)+g0(20,1);
p12_cup = y(t+1,21)+g0(21,1);
p13_cup = y(t+1,22)+g0(22,1);
p14_cup = y(t+1,23)+g0(23,1);
p15_cup = y(t+1,24)+g0(24,1);
p16_cup = y(t+1,25)+g0(25,1);
p17_cup = y(t+1,26)+g0(26,1);
p18_cup = y(t+1,27)+g0(27,1);
p19_cup = y(t+1,28)+g0(28,1);
p20_cup = y(t+1,29)+g0(29,1);
p21_cup = y(t+1,30)+g0(30,1);
p22_cup = y(t+1,31)+g0(31,1);
p23_cup = y(t+1,32)+g0(32,1);
p24_cup = y(t+1,33)+g0(33,1);
p25_cup = y(t+1,34)+g0(34,1);
p26_cup = y(t+1,35)+g0(35,1);
p27_cup = y(t+1,36)+g0(36,1);
p28_cup = y(t+1,37)+g0(37,1);
p29_cup = y(t+1,38)+g0(38,1);
p30_cup = y(t+1,39)+g0(39,1);
p31_cup = y(t+1,40)+g0(40,1);
p32_cup = y(t+1,41)+g0(41,1);
p33_cup = y(t+1,42)+g0(42,1);
p34_cup = y(t+1,43)+g0(43,1);
p35_cup = y(t+1,44)+g0(44,1);
p36_cup = y(t+1,45)+g0(45,1);
p37_cup = y(t+1,46)+g0(46,1);
p38_cup = y(t+1,47)+g0(47,1);
p39_cup = y(t+1,48)+g0(48,1);
p40_cup = y(t+1,49)+g0(49,1);
 
%% Setting the dimension of the matrices 
fx11  = zeros(ny+mx,mx);
fx11p = zeros(ny+mx,mx);
fy   = zeros(mx+ny,ny);
fyp  = zeros(mx+ny,ny);
fx12 = zeros(mx+ny,myx);
 
%% f Function evaluation
% START DISPLAYING f
f(1,1)=exp(-vf_cu)*(exp(d_cu)*(((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^(1 - CHI)/(CHI - 1) - U0d + (PHIzero*exp(n_cu)*(1 - exp(l_cu))^(1 - 1/PHI))/(1/PHI - 1)) - U0 + (AA*BETTA)/exp(evf_cu)^(1/(ALFA - 1))) - 1;
f(2,1)=exp(-w_cu)*(KAPAw*Wss - (Css^CHI0*PHIzero*exp(n_cu)*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^CHI*(KAPAw - 1))/(1 - exp(l_cu))^(1/PHI)) - 1;
f(3,1)=(BETTA*exp(-d_cu)*exp(-pai_cup)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(4,1)=- (exp(-a_cu)*exp(-mc_cu)*exp(w_cu)*exp(l_cu)^THETA)/(Kss^THETA*(THETA - 1)) - 1;
f(5,1)=1 - (exp(-mc_cu)*(ETA + (ZI*exp(pai_cu)*(exp(pai_cu)/PAIss^NU - 1))/PAIss^NU - (BETTA*ZI*exp(-d_cu)*exp(-output_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) - 1))/ETA;
f(6,1)=Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss)) - exp(r_cu);
f(7,1)=(exp(output_cu)*((ZI*(exp(pai_cu)/PAIss^NU - 1)^2)/2 - 1))/(exp(c_cu) + DELTA*Kss) + 1;
f(8,1)=Kss^THETA*exp(-output_cu)*exp(a_cu)*exp(l_cu)^(1 - THETA) - 1;
f(9,1)=exp(-evf_cu)*((exp(vf_cup)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1)))/AA)^(1 - ALFA) - 1;
f(10,1)=exp(-p1_cu)*exp(-r_cu) - 1;
f(11,1)=(BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(12,1)=(BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(13,1)=(BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(14,1)=(BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(15,1)=(BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(16,1)=(BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(17,1)=(BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(18,1)=(BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(19,1)=(BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(20,1)=(BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(21,1)=(BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(22,1)=(BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(23,1)=(BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(24,1)=(BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(25,1)=(BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(26,1)=(BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(27,1)=(BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(28,1)=(BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(29,1)=(BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(30,1)=(BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(31,1)=(BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(32,1)=(BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(33,1)=(BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(34,1)=(BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(35,1)=(BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(36,1)=(BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(37,1)=(BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(38,1)=(BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(39,1)=(BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(40,1)=(BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(41,1)=(BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(42,1)=(BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(43,1)=(BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(44,1)=(BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(45,1)=(BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(46,1)=(BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(47,1)=(BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(48,1)=(BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
f(49,1)=(BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI - 1;
% END DISPLAYING f
 
%% Model derivatives fx11, fx11p, fy, fyp
if getDeriv == 1
% START DISPLAYING fx11
% END DISPLAYING fx11
if t == 1;
    fx11  = zeros(ny+mx,mx);
end
 
% START DISPLAYING fx11p
% END DISPLAYING fx11p
 
 
% START DISPLAYING fy
fy(1,1)=-(exp(-vf_cu)*exp(c_cu)*exp(d_cu))/(Css^CHI0*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^CHI);
fy(2,1)=-(CHI*PHIzero*exp(-w_cu)*exp(c_cu)*exp(n_cu)*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^(CHI - 1)*(KAPAw - 1))/(1 - exp(l_cu))^(1/PHI);
fy(3,1)=(BETTA*CHI*exp(-d_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(5,1)=(exp(-mc_cu)*((B*BETTA*CHI*ZI*exp(-d_cu)*exp(-output_cu)*exp(c_cu)*exp(d_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1)) + (BETTA*CHI*ZI*exp(-d_cu)*exp(-output_cu)*exp(c_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI)))/ETA;
fy(6,1)=PHIc*Rss*exp(c_cu)*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss)) - PHIc_1*Rss*exp(c_ba1p)*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fy(7,1)=-(exp(c_cu)*exp(output_cu)*((ZI*(exp(pai_cu)/PAIss^NU - 1)^2)/2 - 1))/(exp(c_cu) + DELTA*Kss)^2;
fy(11,1)=(BETTA*CHI*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p2_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(12,1)=(BETTA*CHI*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p3_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(13,1)=(BETTA*CHI*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p4_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(14,1)=(BETTA*CHI*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p5_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(15,1)=(BETTA*CHI*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p6_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(16,1)=(BETTA*CHI*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p7_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(17,1)=(BETTA*CHI*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p8_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(18,1)=(BETTA*CHI*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p9_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(19,1)=(BETTA*CHI*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p10_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(20,1)=(BETTA*CHI*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p11_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(21,1)=(BETTA*CHI*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p12_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(22,1)=(BETTA*CHI*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p13_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(23,1)=(BETTA*CHI*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p14_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(24,1)=(BETTA*CHI*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p15_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(25,1)=(BETTA*CHI*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p16_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(26,1)=(BETTA*CHI*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p17_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(27,1)=(BETTA*CHI*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p18_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(28,1)=(BETTA*CHI*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p19_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(29,1)=(BETTA*CHI*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p20_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(30,1)=(BETTA*CHI*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p21_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(31,1)=(BETTA*CHI*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p22_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(32,1)=(BETTA*CHI*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p23_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(33,1)=(BETTA*CHI*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p24_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(34,1)=(BETTA*CHI*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p25_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(35,1)=(BETTA*CHI*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p26_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(36,1)=(BETTA*CHI*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p27_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(37,1)=(BETTA*CHI*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p28_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(38,1)=(BETTA*CHI*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p29_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(39,1)=(BETTA*CHI*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p30_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(40,1)=(BETTA*CHI*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p31_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(41,1)=(BETTA*CHI*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p32_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(42,1)=(BETTA*CHI*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p33_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(43,1)=(BETTA*CHI*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p34_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(44,1)=(BETTA*CHI*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p35_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(45,1)=(BETTA*CHI*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p36_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(46,1)=(BETTA*CHI*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p37_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(47,1)=(BETTA*CHI*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p38_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(48,1)=(BETTA*CHI*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p39_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(49,1)=(BETTA*CHI*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI + (B*BETTA*CHI*exp(-d_cu)*exp(-muz_cup)*exp(-p40_cu)*exp(-pai_cup)*exp(c_cu)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fy(3,2)=(BETTA*exp(-d_cu)*exp(-pai_cup)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(6,2)=-exp(r_cu);
fy(10,2)=-exp(-p1_cu)*exp(-r_cu);
fy(5,3)=-(exp(-mc_cu)*((ZI*exp(2*pai_cu))/PAIss^(2*NU) + (ZI*exp(pai_cu)*(exp(pai_cu)/PAIss^NU - 1))/PAIss^NU))/ETA;
fy(6,3)=PHIpai*Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fy(7,3)=(ZI*exp(output_cu)*exp(pai_cu)*(exp(pai_cu)/PAIss^NU - 1))/(PAIss^NU*(exp(c_cu) + DELTA*Kss));
fy(2,4)=-exp(-w_cu)*(KAPAw*Wss - (Css^CHI0*PHIzero*exp(n_cu)*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^CHI*(KAPAw - 1))/(1 - exp(l_cu))^(1/PHI));
fy(4,4)=-(exp(-a_cu)*exp(-mc_cu)*exp(w_cu)*exp(l_cu)^THETA)/(Kss^THETA*(THETA - 1));
fy(1,5)=(PHIzero*exp(-vf_cu)*exp(d_cu)*exp(l_cu)*exp(n_cu))/(1 - exp(l_cu))^(1/PHI);
fy(2,5)=-(Css^CHI0*PHIzero*exp(-w_cu)*exp(l_cu)*exp(n_cu)*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^CHI*(KAPAw - 1))/(PHI*(1 - exp(l_cu))^(1/PHI + 1));
fy(4,5)=-(THETA*exp(-a_cu)*exp(-mc_cu)*exp(l_cu)*exp(w_cu)*exp(l_cu)^(THETA - 1))/(Kss^THETA*(THETA - 1));
fy(6,5)=PHIl*Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fy(8,5)=-(Kss^THETA*exp(-output_cu)*exp(a_cu)*exp(l_cu)*(THETA - 1))/exp(l_cu)^THETA;
fy(5,6)=-(BETTA*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-output_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(ETA*PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(6,6)=PHIy*Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fy(7,6)=(exp(output_cu)*((ZI*(exp(pai_cu)/PAIss^NU - 1)^2)/2 - 1))/(exp(c_cu) + DELTA*Kss);
fy(8,6)=-Kss^THETA*exp(-output_cu)*exp(a_cu)*exp(l_cu)^(1 - THETA);
fy(4,7)=(exp(-a_cu)*exp(-mc_cu)*exp(w_cu)*exp(l_cu)^THETA)/(Kss^THETA*(THETA - 1));
fy(5,7)=(exp(-mc_cu)*(ETA + (ZI*exp(pai_cu)*(exp(pai_cu)/PAIss^NU - 1))/PAIss^NU - (BETTA*ZI*exp(-d_cu)*exp(-output_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) - 1))/ETA;
fy(1,8)=-(AA*BETTA*exp(-vf_cu)*exp(evf_cu))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*(ALFA - 1));
fy(3,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(5,8)=-(AA*ALFA*BETTA*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-output_cu)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(ETA*PAIss^NU*exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(9,8)=-exp(-evf_cu)*((exp(vf_cup)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1)))/AA)^(1 - ALFA);
fy(11,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(12,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(13,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(14,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(15,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(16,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(17,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(18,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(19,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(20,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(21,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(22,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(23,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(24,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(25,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(26,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(27,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(28,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(29,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(30,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(31,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(32,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(33,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(34,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(35,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(36,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(37,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(38,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(39,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(40,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(41,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(42,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(43,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(44,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(45,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(46,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(47,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(48,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(49,8)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(evf_cu)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1) + 1)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fy(1,9)=-exp(-vf_cu)*(exp(d_cu)*(((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^(1 - CHI)/(CHI - 1) - U0d + (PHIzero*exp(n_cu)*(1 - exp(l_cu))^(1 - 1/PHI))/(1/PHI - 1)) - U0 + (AA*BETTA)/exp(evf_cu)^(1/(ALFA - 1)));
fy(10,10)=-exp(-p1_cu)*exp(-r_cu);
fy(11,11)=-(BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(12,12)=-(BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(13,13)=-(BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(14,14)=-(BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(15,15)=-(BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(16,16)=-(BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(17,17)=-(BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(18,18)=-(BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(19,19)=-(BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(20,20)=-(BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(21,21)=-(BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(22,22)=-(BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(23,23)=-(BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(24,24)=-(BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(25,25)=-(BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(26,26)=-(BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(27,27)=-(BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(28,28)=-(BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(29,29)=-(BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(30,30)=-(BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(31,31)=-(BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(32,32)=-(BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(33,33)=-(BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(34,34)=-(BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(35,35)=-(BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(36,36)=-(BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(37,37)=-(BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(38,38)=-(BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(39,39)=-(BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(40,40)=-(BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(41,41)=-(BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(42,42)=-(BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(43,43)=-(BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(44,44)=-(BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(45,45)=-(BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(46,46)=-(BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(47,47)=-(BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(48,48)=-(BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fy(49,49)=-(BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING fy
 
% START DISPLAYING fx12
fx12(1,1)=(B*exp(-muz_cu)*exp(-vf_cu)*exp(c_ba1)*exp(d_cu))/(Css^CHI0*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^CHI);
fx12(2,1)=(B*CHI*PHIzero*exp(-muz_cu)*exp(-w_cu)*exp(c_ba1)*exp(n_cu)*((exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))/Css^CHI0)^(CHI - 1)*(KAPAw - 1))/(1 - exp(l_cu))^(1/PHI);
fx12(3,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(5,1)=-(B*BETTA*CHI*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-muz_cu)*exp(-output_cu)*exp(c_ba1)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(ETA*PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fx12(6,1)=-PHIc*Rss*exp(c_ba1)*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fx12(11,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(12,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(13,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(14,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(15,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(16,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(17,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(18,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(19,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(20,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(21,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(22,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(23,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(24,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(25,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(26,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(27,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(28,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(29,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(30,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(31,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(32,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(33,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(34,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(35,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(36,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(37,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(38,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(39,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(40,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(41,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(42,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(43,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(44,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(45,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(46,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(47,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(48,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fx12(49,1)=-(B*BETTA*CHI*exp(-d_cu)*exp(-muz_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(c_ba1)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^(CHI - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING fx12
if t < N
% START DISPLAYING fyp
fyp(3,1)=-(BETTA*CHI*exp(-d_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(5,1)=-(BETTA*CHI*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-output_cu)*exp(c_cup)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(ETA*PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1));
fyp(6,1)=PHIc_1*Rss*exp(c_cup)*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fyp(11,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(12,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(13,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(14,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(15,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(16,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(17,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(18,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(19,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(20,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(21,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(22,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(23,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(24,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(25,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(26,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(27,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(28,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(29,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(30,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(31,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(32,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(33,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(34,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(35,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(36,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(37,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(38,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(39,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(40,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(41,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(42,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(43,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(44,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(45,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(46,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(47,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(48,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(49,1)=-(BETTA*CHI*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(c_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^(CHI + 1);
fyp(3,3)=-(BETTA*exp(-d_cu)*exp(-pai_cup)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(5,3)=(exp(-mc_cu)*((BETTA*ZI*exp(-d_cu)*exp(-output_cu)*exp(2*pai_cup)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^(2*NU)*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI) + (BETTA*ZI*exp(-d_cu)*exp(-output_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI)))/ETA;
fyp(6,3)=PHIpai_1*Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fyp(11,3)=-(BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(12,3)=-(BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(13,3)=-(BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(14,3)=-(BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(15,3)=-(BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(16,3)=-(BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(17,3)=-(BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(18,3)=-(BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(19,3)=-(BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(20,3)=-(BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(21,3)=-(BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(22,3)=-(BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(23,3)=-(BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(24,3)=-(BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(25,3)=-(BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(26,3)=-(BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(27,3)=-(BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(28,3)=-(BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(29,3)=-(BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(30,3)=-(BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(31,3)=-(BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(32,3)=-(BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(33,3)=-(BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(34,3)=-(BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(35,3)=-(BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(36,3)=-(BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(37,3)=-(BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(38,3)=-(BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(39,3)=-(BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(40,3)=-(BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(41,3)=-(BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(42,3)=-(BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(43,3)=-(BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(44,3)=-(BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(45,3)=-(BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(46,3)=-(BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(47,3)=-(BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(48,3)=-(BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(49,3)=-(BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(5,6)=(BETTA*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-output_cu)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(ETA*PAIss^NU*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(6,6)=PHIy_1*Rss*exp(PHIl*log(exp(l_cu)/lss) + PHIc*(log(exp(muz_cu)) - log(MUZss) - exp(c_ba1) + exp(c_cu)) + PHIc_1*(log(exp(muz_cup)) - log(MUZss) - exp(c_ba1p) + exp(c_cup)) - PHIpai*(log(exp(paistar_cu)) - log(exp(pai_cu)) + log(PAIss)) - PHIpai_1*(log(exp(paistar_cup)) - log(exp(pai_cup)) + log(PAIss)) + PHIy*log(exp(output_cu)/OUTPUTss) + PHIy_1*log(exp(output_cup)/OUTPUTss));
fyp(3,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(r_cu)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(5,9)=-(AA*ALFA*BETTA*ZI*exp(-d_cu)*exp(-mc_cu)*exp(-output_cu)*exp(-vf_cup)*exp(d_cup)*exp(muz_cup)*exp(output_cup)*exp(pai_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*(exp(pai_cup)/PAIss^NU - 1)*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(ETA*PAIss^NU*exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(9,9)=-(exp(-evf_cu)*exp(vf_cup)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(ALFA - 1))/(AA*((exp(vf_cup)*exp(muz_cup)^((CHI - 1)*(CHI0 - 1)))/AA)^ALFA);
fyp(11,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(12,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(13,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(14,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(15,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(16,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(17,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(18,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(19,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(20,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(21,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(22,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(23,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(24,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(25,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(26,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(27,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(28,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(29,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(30,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(31,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(32,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(33,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(34,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(35,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(36,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(37,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(38,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(39,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(40,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(41,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(42,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(43,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(44,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(45,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(46,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(47,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(48,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(49,9)=-(AA*ALFA*BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(-vf_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^(ALFA - 1))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))*(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI);
fyp(11,10)=(BETTA*exp(-d_cu)*exp(-p2_cu)*exp(-pai_cup)*exp(d_cup)*exp(p1_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(12,11)=(BETTA*exp(-d_cu)*exp(-p3_cu)*exp(-pai_cup)*exp(d_cup)*exp(p2_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(13,12)=(BETTA*exp(-d_cu)*exp(-p4_cu)*exp(-pai_cup)*exp(d_cup)*exp(p3_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(14,13)=(BETTA*exp(-d_cu)*exp(-p5_cu)*exp(-pai_cup)*exp(d_cup)*exp(p4_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(15,14)=(BETTA*exp(-d_cu)*exp(-p6_cu)*exp(-pai_cup)*exp(d_cup)*exp(p5_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(16,15)=(BETTA*exp(-d_cu)*exp(-p7_cu)*exp(-pai_cup)*exp(d_cup)*exp(p6_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(17,16)=(BETTA*exp(-d_cu)*exp(-p8_cu)*exp(-pai_cup)*exp(d_cup)*exp(p7_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(18,17)=(BETTA*exp(-d_cu)*exp(-p9_cu)*exp(-pai_cup)*exp(d_cup)*exp(p8_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(19,18)=(BETTA*exp(-d_cu)*exp(-p10_cu)*exp(-pai_cup)*exp(d_cup)*exp(p9_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(20,19)=(BETTA*exp(-d_cu)*exp(-p11_cu)*exp(-pai_cup)*exp(d_cup)*exp(p10_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(21,20)=(BETTA*exp(-d_cu)*exp(-p12_cu)*exp(-pai_cup)*exp(d_cup)*exp(p11_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(22,21)=(BETTA*exp(-d_cu)*exp(-p13_cu)*exp(-pai_cup)*exp(d_cup)*exp(p12_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(23,22)=(BETTA*exp(-d_cu)*exp(-p14_cu)*exp(-pai_cup)*exp(d_cup)*exp(p13_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(24,23)=(BETTA*exp(-d_cu)*exp(-p15_cu)*exp(-pai_cup)*exp(d_cup)*exp(p14_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(25,24)=(BETTA*exp(-d_cu)*exp(-p16_cu)*exp(-pai_cup)*exp(d_cup)*exp(p15_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(26,25)=(BETTA*exp(-d_cu)*exp(-p17_cu)*exp(-pai_cup)*exp(d_cup)*exp(p16_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(27,26)=(BETTA*exp(-d_cu)*exp(-p18_cu)*exp(-pai_cup)*exp(d_cup)*exp(p17_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(28,27)=(BETTA*exp(-d_cu)*exp(-p19_cu)*exp(-pai_cup)*exp(d_cup)*exp(p18_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(29,28)=(BETTA*exp(-d_cu)*exp(-p20_cu)*exp(-pai_cup)*exp(d_cup)*exp(p19_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(30,29)=(BETTA*exp(-d_cu)*exp(-p21_cu)*exp(-pai_cup)*exp(d_cup)*exp(p20_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(31,30)=(BETTA*exp(-d_cu)*exp(-p22_cu)*exp(-pai_cup)*exp(d_cup)*exp(p21_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(32,31)=(BETTA*exp(-d_cu)*exp(-p23_cu)*exp(-pai_cup)*exp(d_cup)*exp(p22_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(33,32)=(BETTA*exp(-d_cu)*exp(-p24_cu)*exp(-pai_cup)*exp(d_cup)*exp(p23_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(34,33)=(BETTA*exp(-d_cu)*exp(-p25_cu)*exp(-pai_cup)*exp(d_cup)*exp(p24_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(35,34)=(BETTA*exp(-d_cu)*exp(-p26_cu)*exp(-pai_cup)*exp(d_cup)*exp(p25_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(36,35)=(BETTA*exp(-d_cu)*exp(-p27_cu)*exp(-pai_cup)*exp(d_cup)*exp(p26_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(37,36)=(BETTA*exp(-d_cu)*exp(-p28_cu)*exp(-pai_cup)*exp(d_cup)*exp(p27_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(38,37)=(BETTA*exp(-d_cu)*exp(-p29_cu)*exp(-pai_cup)*exp(d_cup)*exp(p28_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(39,38)=(BETTA*exp(-d_cu)*exp(-p30_cu)*exp(-pai_cup)*exp(d_cup)*exp(p29_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(40,39)=(BETTA*exp(-d_cu)*exp(-p31_cu)*exp(-pai_cup)*exp(d_cup)*exp(p30_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(41,40)=(BETTA*exp(-d_cu)*exp(-p32_cu)*exp(-pai_cup)*exp(d_cup)*exp(p31_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(42,41)=(BETTA*exp(-d_cu)*exp(-p33_cu)*exp(-pai_cup)*exp(d_cup)*exp(p32_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(43,42)=(BETTA*exp(-d_cu)*exp(-p34_cu)*exp(-pai_cup)*exp(d_cup)*exp(p33_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(44,43)=(BETTA*exp(-d_cu)*exp(-p35_cu)*exp(-pai_cup)*exp(d_cup)*exp(p34_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(45,44)=(BETTA*exp(-d_cu)*exp(-p36_cu)*exp(-pai_cup)*exp(d_cup)*exp(p35_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(46,45)=(BETTA*exp(-d_cu)*exp(-p37_cu)*exp(-pai_cup)*exp(d_cup)*exp(p36_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(47,46)=(BETTA*exp(-d_cu)*exp(-p38_cu)*exp(-pai_cup)*exp(d_cup)*exp(p37_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(48,47)=(BETTA*exp(-d_cu)*exp(-p39_cu)*exp(-pai_cup)*exp(d_cup)*exp(p38_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
fyp(49,48)=(BETTA*exp(-d_cu)*exp(-p40_cu)*exp(-pai_cup)*exp(d_cup)*exp(p39_cup)*exp(muz_cup)^(CHI*(CHI0 - 1) - CHI0)*(exp(c_cu) - B*exp(-muz_cu)*exp(c_ba1))^CHI*((AA*exp(-vf_cup))/(exp(evf_cu)^(1/(ALFA - 1))*exp(muz_cup)^((CHI - 1)*(CHI0 - 1))))^ALFA)/(exp(c_cup) - B*exp(-muz_cup)*exp(c_cu))^CHI;
% END DISPLAYING fyp
end
end
